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In this work, we investigate the Josephson effect between a two-band superconductor either with the s++ (two 
energy gaps have the same sign and are fully gapped) pairing symmetry or s± (two energy gaps have n phase 
difference and are fully gapped) pairing symmetry and a conventional .s-wave superconductors. The ground 
state, critical current, plasma modes, flux flow dynamics, and response to external ac electric field, possible 
soliton solutions are investigated. For junctions with the charge neutrality breaking, we find a new plasma 
mode for junctions, which gives rise to new resonance peaks in the Josephson flux flow region. Because of the 
frustrated interaction in junctions with s± pairing symmetry, time-reversal-symmetry (TRS) can be broken if 
the frustration is optimized. In the TRS broken (TRSB) state, there is a non-trivial phase difference between 
the two Josephson tunnelling channels, which results in a non-trivial interference. Furthermore, we find a novel 
massless plasma mode at the TRSB transition for junctions with the charge neutrality breaking. In the TRSB 
state, a spontaneous magnetic flux appears where there is a spatial inhomogeneity in the Josephson coupling, 
thus provides a possible smoking-gun evidence for the underlying pairing symmetry. 
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I. INTRODUCTION 

The Josephson effect between two superconductors is a 
hallmark of the macroscopic quantum effect associated with 
superconductivity! 1 ]■ When a phase difference exists between 
two conventional superconductors, spontaneous supercurrent 
flows from one superconductor to the other when they are 
brought together to form a Josephson junction. This is the 
celebrated dc Josephson effect which describes the relation 
between the current and phase difference I s = I c sin(6>i - 82), 
where I c is the critical current and 8, is the superconducting 
phase. When a voltage V is applied to the junction, the super- 
conducting phase rotates according to the ac Josephson effect 
ftd t (8\ - 82) = 2eV. Because of the novel quantum nature of 
the Josephson effect, Josephson junctions have wide applica- 
tions, such as SQUID, electromagnetic devices and quantum 
qubitH). 

The Josephson effect involves the phase of the supercon- 
ductivity therefore depends sensitively on the underlying pair- 
ing symmetry of the superconductors. Thus the Josephson 
effect provides an invaluable tool to pin down the paring 
symmetry of some exotic superconductors, such as Geshken- 
bein, Larkin and Barone phase sensitive measurement for p- 
wave pairing symmetry in heavy -fermion superconductors|3), 
and tri -junction in cuprate superconductors||4). On the other 
hand, the current-phase relation is modified in these exotic 
junctions!?), which broadens our general understanding on 
the Josephson effect and also points a new possible direction 
for applications. 

The multiband superconductivity attracts considerable in- 
terests since the discoveries of MgB 2 superconductor^ and 
iron pnictide superconductors[7]. The pairing symmetry in 
MgB 2 now is well-understood. It has two s-wave energy gaps 
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with the same phase in the cr and n bands, known as s + + 
pairing symmetry! 8 1. While for pnictide superconductors, due 
to vastness of iron pnictide family, discrepancies from differ- 
ent measurements and controversies remain, see Ref.l9 UT31l 
for a review. The so-called extended s-wave or s± pairing 
symmetry with a n phase shift between the hole and electron 
Fermi surface IT61IT71 has attracted considerable attention and 
is favoured by many experiments! 1 8—2 1 ] , and is considered as 
a most probable candidate for iron pnictide superconductors. 

Soon after the discoveries of iron pnictide superconductors, 
Josephson junctions have been fabricated between a conven- 
tional superconductor and these multiband superconductors, 
with an aim to determine the pairing symmetry of pnictide. 
Zhang et. al.\Tl \ fabricated junction between BaKFeAs and 
a conventional superconductor Pb, they clearly observed the 
Fraunhofer pattern in the dependence of the critical current 
on magnetic fields and Shapiro steps when the junction is ir- 
radiated by microwaves. These observations indicate that a 
pure p-wave or of-wave pairing symmetry is unlikely realized 
in pnictide. Similar Fraunhofer pattern is observed by Zhou 
et. al. in a corner junction between BaFeCoAs and Pb[23|. 
Recently, unambiguous half-integer flux quantum jumps have 
been observed in polycrystal NdFesO-Niobium superconduct- 
ing loop|24|, which is a strong evidence for the s± pairing 
symmetry. For a review on MgB 2 junctions and all pnictide 
junctions, please see Ref. ll25H27ll 

From a theoretical point view, for Josephson junctions be- 
tween two-band superconductors and conventional .s-wave 
superconductors, there are two superconducting tunnelling 
channels J s \ sin(#i - 8 S ) and 7 j2 sin(fl 2 - 8 S ), where 8 S is the 
phase of the .s-wave superconductor and 6>,- with i = 1,2 is 
the phase of the two-band superconductor. The relative phase 
between 8\ and 82 is fixed by the underlying pairing symme- 
try. For the s + + pairing symmetry, such as MgB 2 , these two 
channels have the same phase and add constructively, thus the 
phase-current relation is equivalent to a conventional junction 
made of single band superconductors with an effective Joseph- 
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FIG. 1. (color online), (a) Josephson junctions between two-band superconductors with the s + + pairing symmetry and s-wave single 
band superconductors. The interaction between condensates is attractive and they have the same phase in the ground state as shown in (b). 
(c) Josephson junctions between two-band superconductors with the s± pairing symmetry and s-wave single band superconductors. The 
interaction between condensates Oi and 62 is repulsive. Under appropriate conditions, the system is strongly frustrated resulting in TRSB and 
the ground state is depicted in (d). 



son coupling J e g = J sl +J s i- While for the s± pairing symme- 
try, non-trivial phenomena unique to the sign-reversal pairing 
symmetry arise. It was discussed by Agterberg et. al. (28 1 even 
before the discovery of pnictide superconductors that the crit- 
ical current depends non-monotonically on temperatures, and 
may even becomes negative in some temperature region. As- 
suming a n phase shift between the two tunnelling channels, 
numerous phenomena have been demonstrated such as, real- 
ization of the 7rjunction[29, 30], vortex enlargement] 30, new 
Shapiro steps] 32 1, upper bound of the critical current from the 
Ambegaokar - Baratoff relation ll33l . We note that the most 
studies only show quantitatively difference between junctions 
with s + + and s+ pairing symmetry. 

Due to the s+ pairing symmetry, the system is somewhat 
frustrated as shown in Fig. |TJd). Non-trivial phase differ- 
ence between two tunnelling channels beside with an ef- 
fective Josephson coupling J s \ + J& and n with an effective 
Josephson coupling J s \ - is possible. If this happens, 
time-reversal symmetry (TRS) is broken and there are two 
degenerate ground states. The possible violation of TRS is 
first discussed in Ref. [34 1 and later in Ref. lf35l in the context 
of three-band superconductors. As a consequence of time- 
reversal symmetry breaking (TRSB), qualitatively different 
behaviour between junctions made of s+ and s + + supercon- 
ductors may exist, which is hopeful to resolve the dispute of 
the pairing symmetry in pnictide superconductors by observ- 
ing the differences in experiments. 

In this work, we investigate the Josephson effect between 
a two-band superconductor with either s + + or s± pairing 
symmetry and a conventional s-wave superconductor, with an 
emphasis on the possible TRSB. First we derive equations of 
motion for the gauge invariant phase differences based on the 
Lagrangian approach. Then the ground state is obtained by 
minimizing the Josephson energy. We find for junctions with 
s+ pairing symmetry under appropriate conditions, TRSB oc- 
curs. As a consequence of the TRSB, the critical current and 
height of the Shapiro step when junction is shined by an ac 
electric field develop a non-trivial phase dependence on the 
Josephson coupling. We proceed to investigate the Joseph- 



son plasma mode injunctions. For thin superconducting elec- 
trodes, charge neutrality breaking occurs resulting in a mod- 
ified ac Josephson relation. A new plasma dispersion as- 
sociated with an out-of-phase oscillation of phase difference 
emerges. At the TRSB transition, this new plasma mode be- 
comes massless. Because of the existence of the massless 
mode, the Josephson penetration depth diverges and the lower 
critical field H cl of the junction vanishes. In the flux flow re- 
gion, we obtain additional Fiske resonances for a given cavity 
index and an additional Eck resonance because of the charge 
neutrality breaking. 



In the last part of the paper, we discuss possible topological 
excitations in the junction because of the existence of mul- 
tiple degenerate energy minima. In the TRSB state, a new 
type of soliton can be stabilized between the TRSB pair states. 
Finally we discuss appearance of spontaneous magnetic flux 
when there is a spatial variation of the Josephson coupling, 
which occurs only in the junction with s+ pairing symmetry, 
thus provides a possible smoking-gun evidence for the pairing 
symmetry. 



The remaining part of the paper is organized as follows: in 
Sec. II, we derive equations of motion for the gauge-invariant 
phase differences. In Sec. Ill, the boundary condition is de- 
rived. In Sec. IV, we calculate the phase differences at the 
ground state. In Sec. V, the dependence of the critical current 
on magnetic fields is obtained. In Sec. VI, the response of the 
junction to an incident electromagnetic wave is investigated. 
In Sec. VII, we calculate the Josephson plasma modes in the 
junction. In Sec. VIII, the Josephson penetration depth and 
the lower critical field are studied. In Sec. IX, the McCum- 
ber solution is obtained. In Sec. X, we investigate the Fiske 
and Eck resonances in the flux flow region. In Sec. XI, we 
discuss the topological excitations (soliton) in the junction. In 
Sec. XII, possible spontaneous magnetic flux is investigated 
for junctions made of superconductors with s+ pairing sym- 
metry. The paper is closed by conclusions in Sec. XIII. 
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II. MODEL 

We consider a junction between a conventional s-wave su- 
perconductor and a multi-band superconductor either with 
s + + or s± pairing symmetry, as depicted in Fig. [T] The anal- 
ysis can be extended straightforwardly to junctions between 
two two-band superconductors. However the physics is ex- 
pected to be qualitatively the same as the heterotic junctions 
studied here. In the junction between a s + + superconductor 
and a s-wave superconductor, the Josephson couplings among 
different condensates are positive, thus in the ground state they 
have the same phase. While for the junction between a s+ 
superconductor and a .s-wave superconductor, the inter-band 
Josephson coupling between condensates with phase 6\ and 
02 in the s+ superconductor is negative, and the inter-junction 
Josephson coupling between the s-wave superconductor and 
s± superconductor is positive, which results in frustration in 
the system. When the frustration is optimized, different con- 
densates have non-zero phase differences among them which 
indicates the breaking of TRS. We will show below that the 
behavior is qualitatively different for junctions in Fig. [TJa) 
with TRS and Fig. [TJc) with TRSB, which implies a possible 
way to detect the pairing symmetry for multiband supercon- 
ductors. 

The Hamiltonian of the junction is 



< H = <Hi+<H 2 + %, 



(1) 



with the Hamiltonian for the single band superconductor 

tti = f d\ J] 4(r)(e c - nMr) - gcl(r)cl(r)c„(r)c a (r) 

(2) 

and the Hamiltonian for the two-band superconductor 



n 



l,cr 



- 2 f ^r^(r)^_(r)y ;7 ^(r)^(r), (3) 

;,/=i, 2 J 

and the tunnelling between the two superconductors 

h = J] *iA*it + HC > ( 4 ) 

l,CT 

where d, (d^) is the electron creation (annihilation) operator 
in the Z-th band of the two-band superconductor with the dis- 
persion Sdj(k) and the chemical potential /id and spin index 
cr. Vji is the intra-band for / = j and inter-band for I + j scat- 
tering respectively, which can be either repulsive or attractive 
depending, for instance, on the strength of the Coulomb and 
electron-phonon interaction. c] T is the electron creation oper- 
ator in the single band superconductor and g is the electron- 
phonon coupling strength. f/, s is the tunnelling matrix for elec- 
trons between the two superconductors. 

A schematic view of the junction geometry is depicted in 
Fig. [T] External magnetic fields are applied along the y direc- 
tion, which induces inhomogeneous phase differences along 



the x direction. We assume the system is uniform along the 
y direction and the problem becomes two dimensional. We 
proceed to derive equations of motion for the gauge invari- 
ant phase differences for the junction using the Lagrangian 
approach [36, 37 1. The total Lagrangian of system has three 
contributions 



-C = Xi + -Ci + -Lb 



(5) 



with the Lagrangian for the single band superconductor, 
which can be derived from the Hamiltonian Eq. (|2]) using a 
standard method[38 1 



Xi 



d 



A*(r) + ^- c dA(r,t) 
2n 



8t!vI? 



and the Lagrangian the two-band superconductor 



(6) 



i=l,2 



-z 

1=1,2 



A T x (r)--^-dA(r,t) 

Z7T 



+ -z cos(6»i - e 2 ) 

Inc 



(7) 



where fi s and yu,- are the Thomas-Fermi lengths associated with 
charge screening, and A s and Xi are the penetration depths for 
each band, and 9 S and 8, are the superconducting phases for 
different condensates respectively. The effective penetration 
depth for the two-band superconductor is A L 2 = A^ 2 + A^ 2 . Jn 
is the inter-band Josephson coupling, and J\2 < for the s± 
pairing symmetry while J\2 > for the s + + pairing symme- 
try. A* and Aj are electric potentials, and Af and A T X are vec- 
tor potentials at the bottom and top electrodes. (Dq = hc/2e is 
the quantum flux. Here we have introduced the charge energy 
[the first term at the right-hand side of Eqs. (|6]l and ([7])] in the 
superconductors to account for the possible charge neutrality 
brealdng ||3^U4Tl when the thickness of superconducting elec- 
trodes d is comparable to jj s or yu,, which might be realized in 
layered superconductors with strong anisotropy[42|. 
The Lagrangian for the insulating barrier reads 



beb 2 b i 



(8) 



with b being the thickness of the barrier and e c i the dielectric 
constant. The electric field in the barrier is 



I 1 A T — A B 

-d,A b , z - d z Ao = —d t Aba ~ ° r ° - (9) 
c c b 



E b ,z 

and the magnetic field is 

A T - A B 

B b , y = d-A x - d x Ab, z = — ~ ~ d x Ab, z - 
b 

The Josephson coupling Vj is 

Vj = — cos(0 s i) - — cos(0 j2 ), 



2nc 



2nc 



(10) 



(11) 
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with the gauge invariant phase difference 

2nb . 



O 



(12) 



with i = 1, 2 and 6\ - 6 2 — <f> s \ - <f>s2- Here /„■ > are the 
Josephson couplings and can be derived from a microscopic 
theory @ 



J si 



2h 



|A,-A,| _/|A,-[ 



eRbi I A, I + |A,| \|A ; | + |A 



|A S 



(13) 



for a temperature much smaller than the critical temperature, 
where K(x) is the complete eliptic integral of the first kind, 
and A s , A, are the superconducting energy gaps for different 
condensates. The resistance for each channel at the barrier is 



Rbi - 



n 3 



An e 2 Ni(0)N s (0)ti,s' 



(14) 



where Ni(0) is the density of state of quasiparticles in each 
band for the two-band superconductor and N S (Q) is the den- 
sity of state for the single band superconductor. The Joseph- 
son couplings depend on the resistance of the barrier for each 
channel and temperature, which implies a practical way of 
tuning in experiments 1 28 l43ll . 

For convenience, we introduce a length A cl - 
^c<bol{%7T 2 bJ s \) and a frequency cd p1 = c/Cy^/lci). 
We then renormalize the length in unit of A c i and time in 
unit of l/o) p \. Current density is in unit of J s \, electric field 
is in unit of <f>Q0) p i/(2ncb), and magnetic field is in unit of 
® /(2nA cl b). 

We then minimize the Lagrangian in Eq. |5]l using the 
Euler-Lagrange equation. Applying the Euler-Lagrange equa- 
tion with respect to Aq and A B , we have 

J] - [Al + dfif] ~ e h [-bd,A : b - (Aj - A B )\ = 0, (15) 



1=1,2 



— [A B + d,0 s ] + e b [-bd,A\ ~ (4 - A B )\ = 0, (16) 



where a s ^ = u 2 s ^/(db) are parameters characterizing the im- 
portance of the charge neutrality breaking. Equations ( 15 1 and 
( 16 1 can be combined into the following equation using Eqs. 
(9) and ([12) 



±12 a > IA " l " 2 



1 1 

1 

a 1 a 2 



Eb,z = C e Eh, z 



(17) 

Equation ( 17 1 is a modified ac Josephson relation as a result 
of charge neutrality breaking. In the limit a, — > 0, we recover 
the standard ac Josephson relation d,^> S i = £fc,z- We will show 
that a new plasma mode emerges for non-zero a, and a s . 

Similarly minimizing the Lagrangian with respect to A T X 
and A B , we have 

" Z 7 K _ dx0i \ ~ i A * ~ A *- bdxAz ) = ' ( 1 8) 

i=l,2 J 

- - \A B X - d x 6 s ] + (A T X -A B X - bd x A-) = 0, (19) 

where £ sW = A 2 ^ ;) l{db). We then have the relation between 
the magnetic field and spatial derivative of the gauge invariant 



phase difference after rewriting Eqs. ( |T8| > and ( 19 1 using Eqs. 
{TO) and ((T2jl 



Z 

1=1,2 



IT 



1 1 



Bb, y = CbBb,y (20) 



Equation ( 20 1 is a generalization of the phase-magnetic field 
relation in junctions between two single band superconduc- 
tors. In Josephson junctions made of multiband superconduc- 
tors, there exist fractional Josephson vortices according to Eq. 
pO] ), as will be discussed in Sec. XL 

Applying the Euler-Lagrangian equation with respect to 
Ai we obtain the Ampere's law 



d x Bb, y = sin0 s i + 7 s2 sin0 s2 + d,E b ,y 



(21) 



The inter-band Josephson current J 12 sin(^ s i - <p s2 ) does not enter because it does not couple with the gauge field. With the 
help of Eqs. ( 17 1 and ( 20 1, we then arrive at the first equation for the gauge invariant phase difference <p si 



1 



eh 



s " 1 = sin0 sl + / s2 sin0 s2 + — 



(22) 



We still need one more equation for the gauge invariant phase difference. This can be derived by variation of _£ with respect 
to 6j and 6 S , which yields 



— d t Uo(r) + d,G s \ - —d x \A L X - d x s ] + sin0 sI + J s2 sin0 s2 = 0, 

- —8, [A R (r) + d t 6 x \ - -8 X [A R X - d x G x \ - sin^ sI - J l2 sin0 12 = 0, 



- — 8, [A*(r) + d t 6 2 ] - -d x \a r - d x 2 ] - J s2 sin s2 + J l2 sin I2 = 0. 
ebU 2 1 1 (2 1 ' 1 



(23) 
(24) 
(25) 
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Subtracting / £ s x Eq.( 24 1 from Eq. ( 23 1, and using A* (r) + d t 9 s - -a s €bEi„ we obtain another equation for phases 



£s Q.<*1 £, Of! 



Ce \ai a 2 



1 



1 71*- 



i 



+sin0 sl +7 s2 sin0 s2 + — (sin0 s i + J n sin (0 sl - s2 )) = 0. 



(26) 



r 



Equation (21 1 and Eq. (26 1 together with boundary condition 
(derived below) completely describe the phase dynamics in 
the junction. In the limit of no charge neutrality breaking a — > 
0, dt<psi = Ej hz . When \J\2\ » J si , we have (j) s \ = <p S 2 for 
J 12 > 0(s + + pairing symmetry) and <f> s \ = <p S 2 +n for J 12 < 
(i± pairing symmetry). In these limits, the dynamics of the 
junction reduces to a single junction version 



Aid 



J e sin0 s i + 8, 



(27) 



with an effect Josephson coupling J e = 1 + signfVn)/^ and 
an effective penetration depth /l e = 7(£j~* + ( 2 l )/Cb, where 
sign[x] = -1 for x < and sign[x] = 1 for x > 0. Here J e is 
negative for / S 2 < Js\ for s+ pairing symmetry, and Eq. (27 1 



describes the phase dynamics in a n junction[29 30, 43, 44|. 
However, if we express the dynamical equation in terms of 
(p S 2, the sign of the Josephson current becomes positive and 
we have a conventional Josephson junction. 

Dissipations can be introduced through a dissipation func- 
tion 



® = ~B d E 2 b , 



(28) 



where B4 is a damping coefficient. The equation of motion in 
the presence of dissipation can be derived similarly, using the 
Euler-Lagrange equation with dissipation 



69, dx 



6£ 



8{d x 0d 



d_ 

dt 



6D 



(29) 



III. BOUNDARY CONDITION 

Boundary conditions are crucial to determine the dynamics 
inside the junction. For a conventional single band junction, 
the boundary condition is given by d x (p = 27rB a (2A+b)/<$>o be- 
cause the radiation effect is weak|45 1, where B a is the applied 
field, and A is the penetration depth. The boundary condition 
for the multiband junctions cannot be generalized straightfor- 
wardly from single band cases. Using the Usadel equation, 
the boundary condition for multiband superconductors is de- 
veloped in Ref.[46|. In the present work, we assume that 
the supercurrent for each band vanishes at the left and right 
boundary of the two-band superconductor, as shown in Fig. [T| 
and we have 



0>o 
In 



(60 
2n 



d x e 2 =0, 



(30) 



which yields d x <j 



d x (f> S 2 at the boundary using Eq. ( 12 1 



Then from Eq. (20 1, we derive the boundary condition 



d x cf) s i = d x 4> S 2 



<rr 1 + <r 2 



- y (B a ± LJ ext /2) , 



(31) 



where 7 ext is the bias current and L is the length of the junc- 
tion. The second term in the parenthesis at the right hand side 



of Eq. ( 3 1 1 accounts for the magnetic field induced by the bias 
current. The magnetic field inside the junction Bb >y is differ- 
ent from the applied magnetic field B a due to the screening by 
Josephson current. 



IV. GROUND STATE 

In the absence of applied magnetic fields B a — and bias 
current 7 ext = 0, the ground state is determined by minimizing 
the Josephson energy 

E = - cos sl - J s2 cos (pa - J\i cos(0 s i - (f> s2 ), (32) 

which yields 

sin^si + J S 2 sin0 S 2 = 0, (33) 



sin0 s i + J 12 sin(0 s i - <p s2 ) = 0. 



(34) 





\o 2 








FIG. 2. (color online), (a) Phase diagram of Josephson junctions 
with s± pairing symmetry. Here J s2 = J s i = 1. For J n < -0.5, 
TRSB occurs and the transition is continuous. There are two degen- 
erate ground states (j> = ((f> s i,<f>s2) and —f. Only one configuration of 
phases in the ground state with TRSB is shown in the figure, (b) En- 
ergy landscape in the TRSB state with ]\% = -1 and J S 2 = I. There 
are many degenerate and separated energy minima (blue region), 0, 
-<f> and other minima obtained by shifting phases by 2nn with an 
integer n. To contrast the minima, we plot log(£ + 1.6). 
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We see that if = {(p s \ , <f> S 2) is a solution, then 02 = —</> — 
(-<f>si,-<Ps2) is also a solution. If these two solutions are dis- 
tinct, which means that one cannot obtain </>2 from by chang- 
ing the phases by 2nn with an integer n, then TRS is broken 
in the junction. 

In one TRSB state (0ji,0 S 2), the inter-junction Josephson 
current flows from the second band in the two-band supercon- 
ductor to the i-wave superconductor, and then flows back to 
the first band in the two-band superconductor. The net current 
is finally balanced by the interband Josephson current from the 
first band to the second band in the two-band superconductor. 
In the other TRSB state (— tp s \, -tp S 2)> the supercurrent flow re- 
verses direction. The meaning of TRSB becomes explicit in 
this picture. Please note that the supercurrent loop flows in 
the band space, thus is not coupled with the gauge field. No 
spontaneous magnetic flux is induced for spatially homoge- 
neous systems. However as will be demonstrated below, in 
the ground state with TRSB, spontaneous magnetic flux will 
be induced when the Josephson coupling is perturbed locally. 

For J12 > 0, the ground state is trivial (p s \ = <p S 2 = 0. 
However for J\2 < 0, non-trivial solution for S ; occurs due to 
the frustrated interaction. Let us consider J\2 < and I/12I » 
J S i. In this case s j = <p^ + n + 6 with 5« 1. We then expand 



compared to B a . The critical current / is then expressed as 



Eqs. ( 33 1 and ( 34 1 in terms of 6, and obtain 



COS0 s2 = (1 - J S 2)Jn/J S 2, 



(35) 



when J12 > —\Jsz/(l ~ Otherwise <p S 2 — if J s i > 

J s \ and <f> S 2 — n if J S 2 < J s \. Thus TRSB occurs at J 12 = 
— I^s2/(1 - Js2)\- The condition J\2 » is satisfied when J s \ 
and J S 2 are comparable. This is reasonable since frustration is 
maximized when J s \ ~ J S 2- 

For a symmetric coupling = J s \ = 1, we have (p sl = 
—<p S 2 and the ground state can be found exactly. For J 12 > 
-0.5, we have TRS state with <p s \ = <p S 2 = 0. For J\2 < 
-0.5, TRSB occurs with cos <p s \ = -1/(2/12). A typical phase 
diagram and the the energy landscape are shown in Fig. [2] 



/ 1 r L 

- = - (sin0 s i + J s2 sm4> s2 )dx 
L L Jo 
sin® 



[COS 0sio + 7 s2 COS s2 oJ 



(37) 



where we have used the ground state condition Eqs. (33 1 and 
p4] ), and O = ksLll. Besides the conventional factor re- 
sponsible for the Fraunhofer pattern, we have additional term 
accounting for the interference between different tunnelling 
channels. When J\2 > for s + + superconductors, the two 
channels add constructively, while for s± superconductors 
with -J12 » Jsi, we have <p s \ « n + <p S 2, and the two channels 
cancel destructively. The TRSB state with finite < 0„o < x 
interpolates these two limiting cases. We also note that the 
two different TRSB states (f> and -q> has the same critical cur- 
rent. The dependence of critical current on magnetic field is 
depicted in Fig. [3]for three different cases. 



VI. SHAPIRO STEPS 

In this section, we investigate the response of the junction to 
external microwave irradiation. When the external irradiation 
is locked with the internal plasma oscillation, current steps 
are induced known as the Shapiro steps 1 47 1. In the presence 
of incident waves with frequency oj, the voltage across the 
junction can be written as E = Ej c + E sin(wf)- According 
to the generalized ac Josephson relation, the gauge-invariant 
phase differences can be written as 



*si=0sio + <*>ot + A sin(dLit) + 1 



0s2 = 0s2O + o)ot +A sin(wt) + (pi, 



(38) 



(39) 



V. CRITICAL CURRENT 



In this section, we calculate the dependence of the criti- 
cal current on the applied magnetic fields. Because of the 
interference between two tunnelling channels J s \ sin0 s i and 
J S 2 sin0 S 2, the critical current of the junction depends on the 
phase difference between two tunnelling channels. We con- 
sider a short junction L where the screening current can be ne- 
glected. [Precisely, L should be much smaller than the longer 



length in Eqs. (47 1 and (48 1]. The presence of magnetic field 
induces a spatial variation of phases. Without the screening 
effect, the phase can be written as 



+ kgx + (pi\ and . 



>s2 



Ps20 ' 



(36) 



where kg = CbB a / ((1 1 + (2 ') is the phase gradient created 
by the applied field according to the boundary condition Eq. 
(31 1, (pi is the phase created by external current and </> S! -o is 
the phase in the ground state. We have neglected the mag- 
netic field induced by the external current because it is weak 



1 1 1 


1 1 1 1 1 1 1 1 




A J sl +J s^ 




/ \ IJ ^COS^^+J^CQS^^I 




/ \ IJ sl "J s J 




/ 1 V. 1 .V 1 V , A/T 



A5 -10 -5 



FIG. 3. (color online). Dependence of critical current on magnetic 
fields in three different regimes: 7 12 > with (f) sl = (f> s2 , —Jri ^ Jsi 
with <p sl = <p s2 + n and the TRSB state. 
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FIG. 4. (color online). Dispersion of two Josephson plasma modes 
according to Eq. ^45\ and Eq. l |46[ (. We use f = 1, or = 0.05 and 
e b = 1 in the plot. Inset is the energy gap <±>[(k = 0) and w 2 (fc = 0) as 
a function of Jn- The branch (jJi(k) becomes massless at the TRSB 
transition. 



with loq = C e (a^ 1 + a^ 1 ) -1 ^ and A = C e (a^ 1 + a^y 1 E / a>, 
and 4>i an arbitrarily relative phase between the Josephson os- 
cillation and incident wave. Experimentally, when one fixes 
the voltage and tunes the current, 0/ will adjust correspond- 
ingly. The Josephson current then is given by 

/ = <(sin0 s i + J s2 sin0 s2 )>, = 
£j f cos (0sio) + J S 2 cos (0 s2o )} J„(A) sin [(co -nu))t + 7 ] 

n 

(40) 

where J„ is the Bessel function of the first kind and we have 



again used the ground state condition Eqs. (33 1 and (34i 
(■•■), denotes time average. When the resonance condition 



(41) 



C e (a x 1 +a 2 l ) E dc - noj, 

is satisfied, the Shapiro steps appear with the height 

I s = 2 1 [cos (0 slo ) + J s2 cos (0 s2O )] J„(A)\ . (42) 

Again the height of the Shapiro steps depends on the relative 
phase of the different channels, similar to that in the critical 



current in Eq. ( 37 1 



VII. JOSEPHSON PLASMA MODE 

The low energy collective oscillation in a Josephson junc- 
tion is called Josephson plasma, which is composite waves of 
electromagnetic fields and Josephson current. It determines 
lower-energy physical properties and transport properties of 
the junction, such as the Fiske and Eck resonances [ 48 1. More- 
over, the dispersion for the plasma modes can be measured by 
Josephson plasma resonance[49 1, which provides a useful tool 



to extract physical parameters of the junction. In this section, 
we calculate the dispersion of the Josephson plasma in the 
heterotic junction between a multiband superconductor and a 
single band superconductor. 

For convenience of calculations, we take a s — a\ — a 2 — 
(i « 1, (, = (i = (2 - ( and J s \ = J s2 = 1. We expand 
the phases around their ground state value S ;, 0„- = <f> S j + ip S j, 
with j2 = ~<Psi = <Po- We then take the Fourier transform of 
if s i, (fi S i{x, t) ~ (p S i{k, cj) exp[i{kx - cjt)]. Substituting into Eqs. 



( 22 1 and ( 26 1, we have the following equations 



-k 2 a> 2 

+ — - cos 4> Q I (<p sl + ipa) = 0, (43 ) 



-1 , k 2 3f+l 

-0) H — — — — - + 2 cos 4>o + J 12 COS(20()) fsl + 



2e b a £ 3£ + 2 



1 k 2 1 

co 2 - — — — - + cos 0o - J12 cos(20 o ) I f S 2 = (44) 



(2e b a £3£ + 2 
We then derive two branches of the dispersion relation 

2k 2 



(j^[-2 cos 6q + 



3£ + 2' 



+ COS (po + 2J\2 cos(20o) 



(45) 



(46) 



with cos 0o = -l/(27!2). The results are displayed in Fig. [4] 
There are two Swihart velocities with c\ = y/2/(3 + 2g) and 
C2 = -Jeboj^. The splitting of the Josephson plasma into two 
branches thus supports Cherenkov radiation in the junction, 
similar to that in a stack of Josephson junctions [50 51 1. In the 
TRSB state, the plasma modes for the two degenerate TRSB 
pair states and -0 are identical. The mode lo\ corresponds to 
the in-phase oscillation of <p s \ and tp S 2, and is always massive. 
The mode a>2 corresponds to the out-of-phase oscillation, 
which is similar to the Leggett mode ll52l in bulk multiband 
superconductors. In the limit of a — > 0, this mode disappears 
because the out-of-phase oscillation is forbidden according to 
the ac Josephson relation d t <p s \ = d t <p S 2 = Eb tZ . In the presence 
of charge neutrality breaking, the mode 0)2 becomes massless 
at the TRSB transition cos 0o + 2Jn cos(20o) = 0. Close to 
the TRSB transition, the angle between 0,i and S 2 is about 
to grow continuously from to a finite value as shown in Fig. 
|2ja), thus the out-of-phase oscillation becomes a soft mode. 

Near a continuous phase transition, the energy cost for the 
excitation of perturbations is close to zero, thus can be easily 
excited, which in turn tends to destroy the ordered state. At 
the second-order TRSB transition, the out-of-phase oscillation 
s i and S 2 does not cost energy and the phase rigidity for the 
out-of-phase oscillation vanishes. (The phase rigidity for the 
in-phase oscillation is nonzero, see Fig. [4]) In bulk supercon- 
ductors, this phase rigidity is proportional to the superfluid 
density and the mass of photon according to the Anderson- 
Higgs mechanism! 38 1 . The vanish of the phase rigidity in the 
junction means that the photon mass becomes zero. 
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The massless plasma mode at the TRSB transition in Eq. 
( ftp} is remarkable. In a conventional Josephson junction 
where the superconductivity is suppressed in the barrier, the 
mass of the Josephson plasma is reduced compared to that in 
bulk, but the plasma excitation still has an energy gap pro- 
portional to the Josephson coupling. Here we show for the 
first time that there exists genuine massless photon in a junc- 
tion with TRSB. The massless plasma mode is analogous to 
the massless Leggett boson discussed in bulk superconduc- 
tors with TRSB.[53|. The difference is that the Leggett mode 
is a purely phase mode, while the plasma mode is a composite 
mode of superconducting phase and gauge fields. The mass- 
less Josephson plasma mode gives rise novel phenomena as 
will be elucidated in the next section. 



VIII. JOSEPHSON PENETRATION DEPTH AND H c , 

In this section, we calculate the Josephson penetration 
depth and the threshold magnetic field H c \ where magnetic 
flux starts to penetrate into the junctions[54]. The penetration 
depth is determined by the mass of the photon in the junction 
according to the London equation. By setting to — in Eq. 



(45 i and Eq. (46 1, we obtain the Josephson penetration depth 



A n = 1/V(3f + 2)cos( 



A n = 1/ y£ [cos <p + 2J l2 cos(20 o )]- 



(47) 



(48) 



The Josephson penetration depth of the junction is given by 
the larger value of A } \ and Aj2- Being a static property, the 
Josephson penetration depth is independent of a. At the TRSB 
transition, A j2 diverges as a consequence of the vanishing pho- 
ton mass. 

In conventional type II superconductors, magnetic fields 
can penetrate into the superconductor only when they are 
larger than a critical value H c \, due to the screening of fields 
as a consequence of massive photons. This becomes evident 
when we look at H cl ~ <$>o/<l 2 L obtained under the London 
approximation. Here the London penetration depth Ai is in- 
verse proportional to the mass of photon Ai ~ 1/ yfrrip. For 
a conventional Josephson junction, H c \ is reduced as a re- 
sult of weakened superconductivity in the junction, but H c \ 
is nonzero H c \ > 0. In junctions with massless photons at 
TRSB transition, the magnetic field can penetrate easily into 
the system, and the lower critical field becomes zero, H c \ = 0. 



IX. MCCUMBER SOLUTION 

In the absence of applied magnetic fields, the phases are 
uniform along the x direction for a homogeneous junction. 
When the external current exceeds the critical current I c of 
the junction Eq. (36i, the junction switches into a resistive 



state and the phases start to rotate, and such state is called 
the McCumber state [55]. For an overdamped junction (jS^ in 
Eq. (28 i is large » 1. The standard McCumber number 



is given by the junction returns to zero-voltage state 

when the external current is smaller than I c . However, for an 
underdamped junction with j5j « 1, the system remains re- 
sistive, until a current much smaller than I c where the system 
transits into zero-voltage state again[56|. In this section, we 
calculate the IV curve in the McCumber state. 

The phases in the McCumber state, for a symmetric cou- 
pling J S 2 = J si, can be written as 



0s l - 0o + tot + Re [—i<p s \ exp(z'wf)] • 



1>q + <Dt + Re [— i(p S 2 exp(io)t)] . 



(49) 



(50) 



Substituting Eqs. (|49j and ((50]) into Eqs. and (|26j and 
approximating sin(0,i) » — z"exp[z"(0o + tot)], we have 



2e b a 



+ J\2 COS(20o) 



- J 12 COS(20 O ) 



<Ps2 



2e b a 

[2 exp(/0 o ) + exp(-/0 o )] . (51) 



to (fPA + <p s2 ) =4cos(0 o ). 



Equation (51 1 and (52 1 can be solved and we have 



2 cos 0o 2 exp(/0o) + exp(— z'0o) 
sl a) 2 to 2 l(e b a) - 2J\2 cos(20o)' 



</>s2 



2 cos <po 2 exp(/0o) + exp(-z'0o) 



LD- 



u) 2 l{e b a) - 2J 12 cos(20 o ) 



(52) 



(53) 



(54) 



Then the supercurrent contributed from the phase oscillation 
is given by 



-Re 



J s = (sin (0 sl ) + 7 s2 sin (0 s2 )>, 
2 exp(z'^o) + exp(-z'0 o ) 



sin 0o 



w 2 /(e b a) - 2J 12 cos(20 o ) 



(55) 



The supercurrent increases when voltage a> decreases. The to- 
tal current density is J = fya) + J s , with an ohmic contribution 



X. FLUX FLOW REGION: FISKE AND ECK 
RESONANCES 

When strong magnetic fields are applied parallel to the 
junction, magnetic fields penetrate into the junction and form 
Josephson vortices. Under a transport current, the vortices 
are driven by the Lorentz force and move along the junction, 
which induces finite voltage across the junction. Then the 
Josephson plasma is excited. When the plasma is resonant 
with the vortex motion, it induces a large dc part of super- 
current, which manifests as a current step known as the Eck 
step[ 57 ]. For a junction of finite length, the plasma can be res- 
onant with the cavity formed by the junction itself under ap- 
propriate conditions, which also induces current steps known 
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as the Fiske steps[58|. In this section, we calculate the reso- 
nances in the flux flow region. Studies of the flux-flow dynam- 
ics in two-band junctions with an emphasis on the inter-band 
Josephson coupling were presented in Ref. ||59l very recently. 

For a junction with a finite length L and with a symmetric 
coupling J S 2 = J s \, the phase in the flux flow region can be 
written as 

0si — <t>o + ksx + tot + Re [-i<p s \ cos(k m x) exp(/wf)] , (56) 



»s2 - 



• k B x + cot + Re [— i<p S 2 cos(k m x) exp(ia>t)] , (57) 



where kg — Ci,B a /(^ + £T ) is the phase gradient due to the 
applied field B a . Here k m = rrm/L accounts for the geometric 
resonances, and 0o is the ground state value. Substituting Eqs. 
|56]l and ( |57] > into Eqs. ( [22] > and ( po} and using sin(0,i) as 
-/exp[/(0o + k B x + utt)] to the zeroth order, we have 



2e h a 
k 2 



£ 3£ + 2 



J n cos(20 o ) 
1 

<Ps2 = - 



(<Ps2 ~ <fisl) + 



#3f +1 



<fis\ 



(3^ + 2 
[2 exp(z'0 o ) + exp(-/0 o )] F c 



(58) 



+ — (<p sl + <p s2 ) = 2 cos (0 O ) F c , (59) 



F c — — I ofx cos (fc„,x) exp (ik B x) . 
L Jo 



,3£ + 2 2 



where we have introduced the coupling between the flux flow 
and plasma oscillation at the cavity modes k m 



(60) 



(61) 



(62) 



Resonances occur at 



2kl 



3£ + 2 



+ 2/i2 cos(20o), 



that is when the plasma frequency matches the cavity fre- 
quency. Thus in the presence of the charge neutrality break- 
ing, there are two resonant peaks for a given cavity mode m 
because there are two different dispersion branches. The su- 
percurrent induced by the resonance is 



J s = - J (sin (0 s i ) + 7 s2 sin (fisi)), dx 
-iF* 

[ip s i exp(-/0 o ) + J S 2<Ps2 exp(/0 o )] , 



(63) 



where again the terms in the square bracket accounts for the 
interference between two tunnelling channels. Real junction 
inevitably involves dissipation, and the delta-peak resonance 



in Eqs. (61 1 and (62 1 is rounded by the dissipation. A formal 
treatment would introduce dissipation function as discussed in 



Eq. (28 i. Here we introduce phenomenologically the dissipa- 
tion by the replacement u 2 <— to 1 - iB^u). We calculate the IV 




2 4 6 8 10 



FIG. 5. (color online), (a) IV curve in the flux flow region. The 
resonance peaks shown in the figure are caused by the excitation of 
cavity modes (Fiske resonance). Here we only plot the resonance 
mode k„, = n/L with L = 0.2. We use £ = 1, a = 0.05, e b = 1 and 
B a = n. The ohmic contribution in the figure is fjjo) with = 0.01. 
(b) Eck resonance in the flux low region. We use f = 1, a = 0.05, 
e b = 1, B a =7i&ndB d = 0.05. 



curve by solving Eqs. ( 58"1|60"} a nd Eq. (63 1 numerically and 
the result is displayed in Fig. [5|a), where two resonant peaks 
can be clearly observed. 

For a long junction, the geometry resonance is not impor- 
tant. In this case, we need to replace the solutions in Eqs. 
(56 1 and (57 1 by k m <— ks, i.e. the plasma has the same spa- 
tial modulation as the flux flow along the x direction, because 
the plasma oscillation is excited by the flux flow. The expres- 
sions of the resonance conditions and IV curves are the same 
as those of the Fiske resonance except for the replacement of 
Kn <— k B . When the velocity of the flux flow matches the ve- 
locity of plasma (there are two dispersion relations with two 
different velocities), current steps known as Eck steps appear 
as shown in Fig. |5jb). 

For a stack of conventional Josephson junctions, the plasma 
dispersion splits into branches with the number of 
junctions|60-62 1. It gives resonant branches in the TV char- 
acteristic for a given mode. The experimental observation 
of resonant branches thus confirms the plasma splitting in 
junction stacks [63 1. Similarly, the charge neutrality break- 
ing effect discussed here can also be checked experimentally 
by measurement of doubling of resonant branches in the IV 
characteristic. 



XI. TOPOLOGICAL EXCITATIONS 

In the present system, the potential energy has many degen- 
erate minima, which differs in by 27r, as shown in Eq. ( 32 1 



and Fig. |2jb). In the TRSB state, there is additional minima 
at ^ and -0 due to TRSB. Thus the system supports stable 
topological excitations, categorized as solitons or phase kinks 
belonging to the homotopy class of ^o(5' )[64|. There are soli- 
tons with two different types of origins, with one being the 
soliton between the energy minima with 2nn phase shift, and 
the other the soliton between two TRSB pair states, (f> and 
The width of the phase kink is determined by the optimization 
between the potential energy in Eq. ( 32 1 and the energy due 
to the phase gradient. When the width of kink increases, the 
cost in the potential energy increases while the cost in energy 
associated with the phase gradient decreases, and vice versa. 
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The stable kink solution compromises these two energy costs. 
In this section, we investigate the possible soliton solutions in 
the junction. For discussions of phase kink in bulk multiband 
superconductors, please see Ref. ll65H67ll . 

First let us consider the soliton solution by changing (f> s , by 
2n;7T with an integer n,-. The total magnetic flux associated 
with these solitons is given by 



<t(«,) 



(r i l 2n i 7i + ^ l 2n 2 7i)IC b . 



(64) 



It is fractional quantized if n\ + n 2 . Unlike the fractional 
quantized vortices in bulk multiband superconductors ll68ll . the 
energy associated with the fractional soliton in junctions is 
bounded, thus the excitation is thermodynamically stable. We 
calculate numerically a typical configuration of the soliton in 
the junction, where <p s \ changes by 2n and <p s2 does not change 
from the left edge to the right edge. The results are shown in 
Fig. |6|a). The presence of solitons breaks TRS, even though 
in the ground state, the system has TRS. 

Now let us investigate the soliton solution between two 
TRSB pair states (p and -</>, when the ground breaks TRS. 
An exact expression can be found for the symmetric case 
(s = (1 = £2 — £ an d Jsi — J s 2 = 1 ■ In this case 
4>si(x) = -<p s2 {x), and the equation for the spatial variation 
of phase becomes 



d,cf) s i + sin0 s i + J ]2 sin(20 s i) = 0. 



(65) 



We use the Bogomolny method|64| to find the exact solution 



in the following. The energy corresponding to Eq. (65 1 is 

E s = — — ; (cos 4> sl - COS <f> ) - — [cos^O - cos(20 o )] 

(66) 

where cos^o = - 1/(2/12) is the ground state value and we 
introduce them into E s to shift the energy minimum to 0. We 
then consider the following inequality 



1 



V2? 



Vt7 



>o, 



(67) 



with U — ^ V-/12 cos 4>si - 2^=77 ) ' Thus ^ e ener gy of the 



kink has a lower bound 



E, = I- (d,-<Psi f + U> ±^d x <f> sl VC7. (68) 

2 ( y£ 



The lower bound is reached when 



d x <psi = ± y/2£U = + \/2?( V--^i2Cos <f> s \ - — — ) . 

\ 2V--/12/ 

(69) 

Equation (69 1 can be integrated explicitly, which yields a soli- 



ton solution between the two TRSB pair states 
cp s \ = 2tan~' 



1 - 2J 12 



tanh 



•47? 



2 V^V 3 ^ 



(70) 



where + corresponds to soliton and - corresponds to anti- 
soliton. The energy associated with the kink is 



E s 



12 sin 0o 



1 



= 00 



(71) 



In the symmetric case as we studied here, there is no magnetic 
flux associated with the soliton because the phase gradient in 
each channel compensates exactly d x (<f> s \ +4> s i) = 0. However, 
for more general cases, soliton carries finite magnetic flux. We 
calculate numerically the soliton solution and the associated 
magnetic flux for asymmetric parameters £\+ + ( s , and the 
results are displayed in Fig. |6|b). Suppose at the left domain, 
the phases belong to one ground state -^10 and — <f> S 2Q, and at 
the right domain, the phases belong to the other ground state 
0,io and <p S 2o, the total magnetic flux is quantized in terms of 



2(£, V.sl0 + £ 2 Vs20)/Ci. 



(72) 



The discussions of the solitons so far are restricted to the 
static case where the stability of solitons are guaranteed by 
the energy landscape of the Josephson coupling. In the dy- 
namic region when a transport current is present, solitons can 
be stabilized with some voltages. Solitons can also be created 
due to the dynamic instability of superconducting phase in the 
Josephson junctions. |69] The solitons start to move because 
of the Lorentz force. When the velocity of soliton matches 
the velocity of plasma, current steps known as the zero-field 
steps arise at these voltages [70]. Due to the existence of vari- 
ous types of fractionally quantized solitons and the existence 
of two Swihart velocities for plasma in the present system, we 
expect appearance of many zero-field steps in the IV curve. 

The existence of solitons associated with flux in Eq. ( 64 1 



does not depends on the pairing symmetry of the two-band 
superconductors, while the soliton between two TRSB pair 
states in Eq. f70| ) exists only in superconductors with s+ pair- 
ing symmetry, where TRSB in the ground state is possible. 



XII. SPONTANEOUS MAGNETIC FLUX 

In this section, we investigate the spontaneous magnetic 
flux in junctions unique to superconducting electrodes with 
the s+ pairing symmetry, which points a practically useful 
way to detect the pairing symmetry. We consider a junction 
with spatial variation of the Josephson couplings, which is 
likely realized in real systems due to inhomogeneity, such as 
variation of the thickness of barrier, or can be tuned intention- 
ally as in Ref. l24l . Let us first consider a step modulation 
of the Josephson coupling, which is amenable to analytical 
calculation 



J S 2 ~ Jsi = J'[2®(X) - 1], 



(73) 



with the Heaviside step function ©(x). We further assume 
-J12 s> J si , and the spatial variation of phase is given by 



according to Eq. (27 1 



Xfifa - J'[-2®(x) + 1] sin0 sl = 0, 



(74) 
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FIG. 6. (color online), (a) Profile of phases and magnetic field in 
(a) a conventional 2n solition, where <f> sl changes by 2n, (b) soliton 
between two TRSB pair states. Here f, = 10, <f 2 = 24, & = 20, 
J,\ = 1.0, / s2 = 1.5. J u = 1.0 in (a) and J n = -1.0 in (b). 



with the boundary condition d x <p s \{x = +L/2) = 0. The length 
of the junction L is assumed to be small L <K A e . In this case, 
the solution can be expressed as 



(75) 



<f>sl(x) = </) s io + T](x)shl(f) sW , 



with rj{x) « 1. Substituting the solution into Eq. ( |74| i, we 
obtain the equation for rj(x) 

A 2 e d 2 x Tj(x) - J'(-2@(x) + 1) [1 + rj(x) cos jlo ] = 0. (76) 

The first term in the square bracket dominates and to the first 
order, we can safely neglect the second term in the square 
bracket. Integrating along the x direction yields the magnetic 




FIG. 7. (color online). Profile of the phases and magnetic field with 
a step modulation of the losephson coupling in a junction with s± 
pairing symmetry. Here 0,2 = <psi + tc and X e = J' = 1. 



flux 



B 



C b 



sii\(f> s \od x T](x) = /'sin( 



(77) 



cpsio is given by the condition that spatial average of Eq. (76 1 
vanishes because of the boundary condition d x <p s \ — as no 
external magnetic field is applied 



r. 

J-L/2 



J'[-2®(x) + 1]tj(x) cos ip s iodx = 0, 



(78) 



which yields s io = n/2. 

The magnetic flux inside depends on the length of junc- 
tion and is not quantized for a short junction L <K A e . How- 
ever for a long junction L » A e , magnetic flux only occurs 
in the region where the Josephson coupling changes and it is 
quantized due to the complete screening. We investigate the 
magnetic flux in a long junction L » A e , where the nonlinear 
effect of Josephson current must be treated self-consistently. 
The solution can be constructed as follows: for x < 0, 



: 4 arc tan 



exp 



( V7 7 



(x + x ) 



and for x > 



A* = 4 arctan 



exp 



— (x- Xq) 

A e 



(79) 



(80) 



At x = 0, the derivative of <p sl is automatically continu- 
ous under this construction, xo is determined by the con- 
tinuity condition <p^(x = 0) = <p^{x = 0), which yields 

Xo = /l e ln(— 1 + V2) / yfl' . The total magnetic flux inside 
the junction is given by 



rL/2 
J-L/2 



Bdx — 



c b 



I d x <f> s i = 
Jo 



c b 



n, (81) 



is quantized since <f> s \ run from at x — -L/2 to n at x — L/2. 
The profile of phases and magnetic field are shown in Fig. [7] 
The modulation of the Josephson coupling Eq. ( 73 1 can be 
simplified into a setup shown in Fig. [8|a), where a loop made 
of conventional s-wave superconductor is in contact with a 
two-band superconductor and forms two point junctions. This 
geometry is the same as the one studied experimentally in Ref. 
||29ll . When \J 12 \ » /.„, Q\ = n + 6 2 . If at the left contact 
■/si > Js2> tnen = Q\\ However, if one can achieve at the 
right end < J s i, then 6f — 02 = 6\ + n when the self- 
inductance is small for a large loop. Then there is a ji phase 
shift between the two ends of the s-wave superconductor. This 
is a new way to realize the n junctional | by inserting a two- 
band superconductor with s± pairing symmetry between two 
conventional s-wave superconducting electrodes. For a super- 
conducting loop containing a n junction, TRSB occurs with 
a spontaneous magnetic flux of <&o/2. If this spontaneous 
magnetic flux is measured, then it can unequivocally prove 
the sign-reversal pairing symmetry for the two-band super- 
conductors. 
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FIG. 8. (color online), (a) Experimental proposal to detect the 
pairing symmetry in multiband superconductors, (b) Grain boundary 
Josephson junctions in multiband superconductors. 



The above analysis can be generalized into Josephson junc- 
tions in grain boundaries of s+ multiband superconductors. 
As shown in Fig. [HJb), there are four tunnelling channels in 
the grain boundary junction, with two diagonal and two off- 
diagonal tunnellings. The total energy can be written as 



cos(0,i - 

-7,21 COS(0,2 - 

-7, cos(0,i 



m)- 
m)- 
0a) 



7,i2 cos(6»,i 

7,22 008(0,2 ' 
- Jb COS(0 n - 



e b2 ) 
e b2 ) 



(82) 



where J la p > with a, ft = 1 , 2 are the inter-junction cou- 
plings, and 7, < 0, 7i < are inter-band couplings. For |7,| 
and \J b \ » /fojg > 0, 6,\ — t2 + it and 6 b \ = &b2 + x, then the 
Josephson energy is simplified into 



E = (-7 f n - 7,22 + 7,12 + 7 f2 1) COS(6»,2 - 6 b2 ). 



(83) 



If at some part of the junction (-7,n - 7,22 + 7,i2 + 7,2i) > 0, 
while at other part (— 7,i i - 7,22 + Jai + Ja\) < 0, then we have 
exactly the same situation as studied in Eqs. f73j ) and f74} . 
Spontaneous magnetic flux will appear in the grain boundary 
if this condition holds. 

We then investigate the response of a junction in the TRSB 
state when spatial inhomogeneities of the Josephson cou- 
plings are present. In the TRSB state, S( - depends on the 
Josephson couplings. The inhomogeneities induce variation 
of <f> S j, and thus create magnetic field. We assume a weak de- 
fect in the 7 S 2 channel 7 S 2 = 7 S 2 + A S 2(x), and the response is 
<f>si = 0si + (psi wim A S 2 <s 1 and <p S i <K 1. For simplicity, we 
assume two channels are identical £i = £z = £ s = 7;2 = J s i 
and S 2 — -(f> s \ — <po- Substituting these expressions into Eqs. 



1 1 1 
"J^l.O-O.lSCx) 


— 








FIG. 9. (color online). Profile of the phases and magnetic field 
with a defect in the junction with s± pairing symmetry obtained from 
Eqs. ^88|90| >. The defect is modelled as inhomogeneous Josephson 
coupling in the junction, (a) and (b) correspond to two TRSB pair 
states. Inset is the phase configuration in the ground state. 



( 22 1 and ( 26 1 and then expanding to the linear order, we have 



3<T + 2 



cos 4> (p si + A s2 (x) sin O + J s2 tp s2 cos </> , 

(84) 

34" + 1 d 2 x ifi s i dl<p s2 . 

1 h 2 cos d>o<p s \ + A s2 (x) sin <bn 

3£ + 2 ( £(3^ + 2) v W W 

+7 s2 COS O ^s2 + J 12 COS(20 O ) (^sl - Vs2) = 0. 

(85) 



Equations ( 84 1 and ( |85] l can be solved in the Fourier space and 
the solution is given by 

sin O \2k\\ + + P 2 (2 + 30 - q 2 d 

= WTp^¥) m > m 

sin 0o \2k 2 (l + 20 + p 2 (2 + 30 + q 2 (] 
2 { k 2 +P 2 ){ k 2 + q 2 ) m ' (87) 

with p 2 = [cos 0o + 2 cos (20o) Jn] £, and q 2 — cos 0o(2 + 3£). 
We assume a point defect A(x) = A^6(x). Then we obtain the 
induced modulation of phase in real space 

P -\Ap n 



<p s l(x) = - 



Ao sin < 



e -qlxl q 

cos 0o cos 00+2 cos (20o) J 12 



(88) 




FIG. 10. (color online). Profile of the phases and magnetic field 
with a Gaussian-type defect Ja(x) = 1 - 0.4exp(-(x - L/2) 2 ) in the 
Josephson coupling in a junction with s± pairing symmetry obtained 
by numerical calculations, (a) and (b) correspond to two TRSB pair 
states. Inset is the phase configuration in the ground state. For clarity, 
the magnetic field is amplified by 20 times. 
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<Psl(x) = ~ 



A o sin0o 



q 



e -WPp 



cos 0o cos <po + 2 cos (20o ) ^12 



and the associated magnetic flux is given by 



B 



2Ao sin0osign[jc] e ' A ' 9 ^ 2 



3^ + 2 



4 cos ( 



(89) 



(90) 



Please note that the magnetic flux is singular at x = because 
of the 6 function used for the defects. Since Eqs. ( [88H90| > 
all contain sin 0o> the response to the same defect for distinct 
TRSB state differs by a sign, as shown in Fig. [9] (a) and (b). 
This can be understood by looking at the phases of the three 
condensates. We have the freedom to fix the phase of the s- 
wave superconductor 8 S = 0. As shown in the inset of Fig. 
|9|a), when is suppressed due to the local defect, the attrac- 
tion between 8 S and 82 decreases, thus 82 decreases (rotates 
leftwards). As a result 8\ also decreases because the repul- 
sion between 8\ and 82 is reduced. The width of the region of 
phase variation is optimized by energy cost due to the phase 
gradient. While for the other TRSB ground state shown in Fig. 
[9jb), when decreases, 82 increases (rotates rightwards) as a 
result of the reduced attraction. Meanwhile 8\ also increases. 

For a strong defect, we calculate the phases and magnetic 
field numerically. We model the point defect by a Gaussian 
distribution, and the results are presented in Fig. 10 Mag- 



netic flux is induced near the defect but the integrated mag- 
netic flux in the junction vanishes. Meanwhile the response 
for two distinct TRSB states is different, consistent with the 
analytical results in Eqs. ([88](90jl. 



XIII. CONCLUSIONS 

In this work, we have studied systematically the Josephson 
effect between a two-band superconductor and a conventional 
i-wave superconductor. We consider both the s + + and s+ 
pairing symmetries for the two-band superconductor. Due to 
the multiband nature of the superconducting electrode, there 
are two tunnelling channels in the junction, which gives rise to 
complicated interference in physical quantities, depending on 
the underlying pairing symmetry. Moreover for junctions with 
s+ pairing symmetry, there exists frustrated interaction among 



different condensates and under appropriate conditions, time- 
reversal symmetry is broken. Depending on the competition 
of inter-junction Josephson couplings and inter-band Joseph- 
son coupling, the interference between the two tunnelling 
channels can change continuously from adding constructively 
where they have the same phase to cancelling destructively 
where they have a n phase shift. The interference manifests 
itself in the critical current, Shapiro steps, Fiske current steps 
and Eck current steps. 

In the case of thin superconducting electrodes, charge neu- 
trality can be broken in the electrodes. Because of the charge 
neutrality breaking, out-of-phase oscillations of the gauge- 
invariant phase difference between two channels are possible, 
which gives a new plasma mode, in addition to the in-phase 
plasma oscillations. The energy gap of the out-of-phase mode 
vanishes at the time-reversal symmetry breaking transition, 
which is a first example of massless plasma mode in super- 
conductors. Because of the existence of two plasma modes, 
there are additional Fiske and Eck resonances in the Joseph- 
son flux flow region. 

A long junction supports topological excitation of solitons. 
When the gauge invariant phase difference associated with 
one tunnelling channel changes by multiple 2n that is dif- 
ferent from the change of phase difference in the other chan- 
nel, fractional quantized magnetic solitons are stabilized in the 
junction. In contrast to the fractional vortices in bulk multi- 
band superconductors, the fractional solitons have finite en- 
ergy thus are thermodynamically stable. For junctions with 
time-reversal symmetry breaking, a new type of soliton exci- 
tations can be created between two distinct time-reversal pair 
states. 

Finally we showed that injunctions with time-reversal sym- 
metry breaking, disorders induce magnetic flux in the junc- 
tions, which points a unique way to detect the pairing symme- 
try of multiband superconductors. 



XIV. ACKNOWLEDGEMENT 

The author is grateful to L. N. Bulaevskii for helpful dis- 
cussions. This work was supported by the US Department of 
Energy, Office of Basic Energy Sciences, Division of Materi- 
als Sciences and Engineering. 



[1] B. D. Josephson, Phys. Lett. 1, 251 (1962). 

[2] A. Barone and G. Paterno, Physics and Applications of The 

Josephson Effect (Wiley, 1982). 
[3] V. B. Geshkenbein, A. I. Larkin, and A. Barone, Phys. Rev. B 

36, 235 (1987). 

[4] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 (2000). 
[5] A. A. Golubov, M. Y. Kupriyanov, and E. Il'ichev, Rev. Mod. 

Phys. 76, 411 (2004). 
[6] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and 

J. Akimitsu, Nature 410, 63 (2001). 
[7] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. 

Chem. Soc. 130, 3296 (2008). 



[8] X. X. Xi, Rep. Prog. Phys. 71, 116501 (2008). 
[9] K. Ishida, Y. Nakai, and H. Hosono, J. Phys. Soc. Jpn 78, 
062001 (2009). 
[10] J. Paglione and R. L. Greene, Nat. Phys. 6, 645 (2010). 
[11] D. C. Johnston, Adv. Phys. 59, 803 (2010). 
[12] H.-H. Wen and S. Li, Annu. Rev. Condens. Matter Phys. 2, 121 
(2011). 

[13] F. Wang and D.-H. Lee, Science 332, 200 (201 1). 

[14] P. J. Hirschfeld, M. M. Korshunov, and 1. 1. Mazin, Rep. Prog. 

Phys. 74, 124508 (2011). 
[15] G. R. Stewart, Rev. Mod. Phys. 83, 1589 (2011). 



14 



[16] 1. 1. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, Phys. 

Rev. Lett. 101, 057003 (2008). 
[17] K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. Kontani, 

and H. Aoki, Phys. Rev. Lett. 101, 087004 (2008). 
[18] H. Ding, P. Richard, K. Nakayama, K. Sugawara, T. Arakane, 

Y. Sekiba, A. Takayama, S. Souma, T. Sato, T. Takahashi, 

Z. Wang, X. Dai, Z. Fang, G. F. Chen, J. L. Luo, and N. L. 

Wang, Europhys. Lett. 83, 47001 (2008). 
[19] L. Wray, D. Qian, D. Hsieh, Y. Xia, L. Li, J. G. Checkelsky, 

A. Pasupathy, K. K. Gomes, C. V. Parker, A. V. Fedorov, G. F. 

Chen, J. L. Luo, A. Yazdani, N. P. Ong, N. L. Wang, and M. Z. 

Hasan, Phys. Rev. B 78, 184508 (2008). 
[20] K. Terashima, Y. Sekiba, J. H. Bowen, K. Nakayama, T. Kawa- 

hara, T. Sato, P. Richard, Y.-M. Xu, L. J. Li, G. H. Cao, Z.-A. 

Xu, H. Ding, and T. Takahashi, Proc. Natl. Acad. Sci. USA 

106, 7330 (2009). 
[21] Y. Zhang, L. X. Yang, F. Chen, B. Zhou, X. F. Wang, X. H. 

Chen, M. Arita, K. Shimada, H. Namatame, M. Taniguchi, J. P. 

Hu, B. P. Xie, and D. L. Feng, Phys. Rev. Lett. 105, 117003 

(2010). 

[22] X. H. Zhang, Y. S. Oh, Y. Liu, L. Q. Yan, K. H. Kim, R. L. 

Greene, and I. Takeuchi, Phys. Rev. Lett. 102, 147002 (2009). 
[23] Y. R. Zhou, L. Y. R., J. W. Zuo, R. Y. Liu, S. K. Su, G. F. Chen, 

J. L. Lu, N. L. Wang, and Y. P. Wang, arXiv:0812.3295 (2008). 
[24] C. T. Chen, C. C. Tsuei, M. B. Ketchen, Z. A. Ren, and Z. X. 

Zhao, Nature Phys. 6, 260 (2010). 
[25] A. Brinkman, D. Mijatovic, H. Hilgenkamp, G. Rijnders, 

I. Oomen, D. Veldhuis, F. Roesthuis, H. Rogalla, and D. H. A. 

Blank, Supercond. Sci. Technol. 16, 246 (2003). 
[26] X. X. Xi, Supercond. Sci. Technol. 22, 043001 (2009). 
[27] P. Seidel, Supercond. Sci. Technol. 24, 043001 (2011). 
[28] D. F. Agterberg, E. Dernier, and B. Janko, Phys. Rev. B 66, 

214507 (2002). 

[29] W. Q. Chen, F. J. Ma, Z. Y. Lu, and F. C. Zhang, Phys. Rev. 

Lett. 103, 207001 (2009). 
[30] D. Parker and 1. 1. Mazin, Phys. Rev. Lett. 102, 227007 (2009). 
[31] Y. Ota, M. Machida, T. Koyama, and H. Matsumoto, Phys. Rev. 

B 81, 014502 (2010). 
[32] Y. Ota, M. Machida, and T. Koyama, Phys. Rev. B 82, 

140509(R) (2010). 
[33] Y. Ota, N. Nakai, H. Nakamura, M. Machida, D. Inotani, 

Y. Ohashi, T. Koyama, and H. Matsumoto, Phys. Rev. B 81, 

214511 (2010). 

[34] T. K. Ng and N. Nagaosa, Europhys. Lett. 87, 17003 (2009). 

[35] V. Stanev and Z. Tesanovic, Phys. Rev. B 81, 134522 (2010). 

[36] E. Simanek, Inhomogeneous Superconductors: Granular and 
Quantum Effects (Oxford University Press, New York, 1994). 

[37] Y. Ota, M. Machida, T. Koyama, and H. Matsumoto, Phys. Rev. 
Lett. 102, 237003 (2009). 

[38] A. Alexander and B. D. Simons, Condensed Matter Field The- 
ory (Cambridge University Press, Cambridge, 2010). 

[39] H. Suhl, Phys. Rev. Lett. 14, 226 (1965). 

[40] T. Koyama and M. Tachiki, Phys. Rev. B 54, 16183 (1996). 

[41] M. Machida, T. Koyama, and M. Tachiki, Phys. Rev. Lett. 83, 
4618 (1999). 



[42] R. Kleiner, F. Steinmeyer, G. Kunkel, and P. Miiller, Phys. Rev. 

Lett. 68, 2394 (1992). 
[43] J. Linder, I. B. Sperstad, and A. Sudbo, Phys. Rev. B 80, 

020503 (2009). 

[44] W.-Q. Chen and F.-C. Zhang, Phys. Rev. B 83, 212501 (201 1). 
[45] L. N. Bulaevskii and A. E. Koshelev, Phys. Rev. Lett. 97, 
267001 (2006). 

[46] A. Brinkman, A. A. Golubov, and M. Y. Kupriyanov, Phys. 

Rev. B 69, 214407 (2004). 
[47] S. Shapiro, Phys. Rev. Lett. 11, 80 (1963). 
[48] M. Cirillo, N. Gronbech- Jensen, M. R. Samuelsen, M. Salerno, 

and G. V. Rinati, Phys. Rev. B 58, 12377 (1998). 
[49] A. J. Dahm, A. Denenstein, T. F. Finnegan, D. N. Langenberg, 

and D. J. Scalapino, Phys. Rev. Lett. 20, 859 (1968). 
[50] R. G. Mints and I. B. Snapiro, Phys. Rev. B 52, 9691 (1995). 
[51] G. Hechtfischer, R. Kleiner, A. V. Ustinov, and P. Miiller, Phys. 

Rev. Lett. 79, 1365 (1997). 
[52] A. J. Leggett, Prog. Theor. Phys. 36, 901 (1966). 
[53] S. Z. Lin and X. Hu, Phys. Rev. Lett. 108, 177005 (2012). 
[54] M. Tinkham, Introduction to Superconductivity (McGraw-Hill, 

Inc., New York, 1996). 
[55] R. Kleiner, T. Gaber, and G. Hechtfischer, Phys. Rev. B 62, 

4086 (2000). 

[56] X. Hu and S. Z. Lin, Supercond. Sci. Technol. 23, 053001 
(2010). 

[57] R. E. Eck, D. J. Scalapin, and B. N. Taylor, Phys. Rev. Lett. 13, 
15 (1964). 

[58] M. D. Fiske, Rev. Mod. Phys. 36, 221 (1964). 

[59] J. H. Kim, B.-R. Ghimire, and H.-Y. Tsai, Phys. Rev. B 85, 

134511 (2012). 
[60] K. L. Ngai, Phys. Rev. 182, 555 (1969). 

[61] S. Sakai, P. Bodin, and N. F. Pedersen, J. Appl. Phys. 73, 2411 
(1993). 

[62] S. Sakai, A. V. Ustinov, H. Kohlstedt, A. Petraglia, and N. F. 

Pedersen, Phys. Rev. B 50, 12905 (1994). 
[63] A. V. Ustinov, H. Kohlstedt, M. Cirillo, N. F. Pedersen, G. Hall- 

manns, and C. Heiden, Phys. Rev. B 48, 10614 (1993). 
[64] N. Manton and P. M. Sutcliffe, Topological solitons (Cambridge 

University Press, Cambridge, 2004). 
[65] Y. Tanaka, Phys. Rev. Let. 88, 017002 (2002). 
[66] S. Z. Lin and X. Hu, New J. Phys. 14, 063021 (2012). 
[67] V. Vakaryuk, V. Stanev, W. Lee, and A. Levchenko, 

arXiv: 1203.4554 (2012). 
[68] E. Babaev, Phys. Rev. Lett. 89, 067001 (2002). 
[69] S. Pagano, M. P. Soerensen, R. D. Parmentier, P. L. Chris- 
tiansen, O. Skovgaard, J. Mygind, N. F. Pedersen, and M. R. 

Samuelsen, Phys. Rev. B 33, 174 (1986). 
[70] T. A. Fulton and R. C. Dynes, Solid State Commun. 12, 57 

(1973). 

[71] L. N. Bulaevskii, V. V. Kuzii, and A. A. Sobyanin, JETP Lett. 
25, 290 (1977). 



